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Abstract: We perform gfobal fits to the parameters of tlie Constrained Minimal Super- 
symmetric Standard Model (CMSSM) and to a variant with non-universal Higgs masses 
(NUHMl). In addition to constraints from low-energy precision observables and the cos- 
mological dark matter density, we take into account the LHC exclusions from searches in 
jets plus missing transverse energy signatures with about 5fb~^ of integrated luminosity. 
We also include the most recent upper bound on the branching ratio Bg — )• from LHCb. 
Furthermore, constraints from and implications for direct and indirect dark matter searches 
are discussed. The best fit of the CMSSM prefers a light Higgs boson just above the ex- 
perimentally excluded mass. We find that the description of the low-energy observables, 
{g — 2)^ in particular, and the non-observation of SUSY at the LHC become more and 
more incompatible within the CMSSM. A potential SM-like Higgs boson with mass around 
126 GeV can barely be accommodated. Values for 13{Bs — )■ nn) just around the Standard 
Model prediction are naturally expected in the best fit region. The most-preferred region 
is not yet affected by limits on direct WIMP searches, but the next generation of experi- 
ments will probe this region. Finally, we discuss implications from fine-tuning for the best 
fit regions. 
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1 Introduction 

Supersymmetry (SUSY) [1] is one of the theoretically best-motivated extensions of the 
Standard Model of particle physics (SM). Furthermore, there are intriguing complementary 
experimental findings: the measured coupling constants of the electro-magnetic, weak, and 
strong forces unify when extrapolated to high energies using the renormalization group 
equations (RGEs) of SUSY, but not with the RGEs of the SM [2]. The observed dark 
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matter (DM) density in the Universe [3] can easily be explained by a stable neutralino 
lightest supersymmetric particle (LSP), acting as a weakly interacting massive particle 
(WIMP) [4]. The measured anomalous magnetic moment of the muon [5] appears to 
deviate slightly from its SM prediction [6] but can be naturally accommodated within 
SUSY with relatively light SUSY particles. Moreover, several precision observables like the 
mass of the W boson and the effective weak mixing angle sin^ 6^^, and some rare B-meson 
decay branching ratios are also sensitive to the parameters of SUSY models. On the other 
hand, SUSY particles have not been directly observed at colliders [7-38] and there is no 
convincing evidence for the existence of WIMPs from direct or indirect searches [39-42] . 

The simplest supersymmetric extension of the SM is the Minimal Supersymmetric 
Standard Model (MSSM) [43-45]. Here the proton is stable due to conserved R-parity, 
which also guarantees the stability of the LSP. In its most general form (but with massless 
neutrinos) the MSSM has 124 free parameters, including the 19 of the SM.^ However, 
current data is insufficient to place meaningful constraints on this general model. For 
phenomenological studies it is thus essential to consider more restrictive models. The most 
widely considered is the constrained MSSM (CMSSM) [43, 45] which has only 5 new free 
parameters beyond the SM^ 

Mo, Ml/2, ^0, tan/3, sgn(^). (1-1) 

Here Mq and Mi ^2 denote the universal soft supersymmetry breaking scalar and gaugino 
masses at the unification scale, respectively. ^0 is the universal soft supersymmetry break- 
ing trilinear scalar coupling, tan /3 is the ratio of the vacuum expectation values of the two 
CP-even neutral Higgs fields and n is the Higgs mixing parameter in the superpotential. 

Ignoring the direct searches at the ATLAS and CMS experiments at the Large Hadron 
Collider (LHC) for now, all data is compatible with the 5 parameters of this highly simpli- 
fied model (see e.g. [51-56]). The strictest constraints come from interpreting the WMAP 
measurement in terms of the neutralino relic density [3] and from the anomalous magnetic 
moment of the muon. The best fits prefer Mq ^ 100 GeV and M1/2 ^ 400 GeV [51-56], see 
also [57]. 

In a minimal non-universal Higgs mass model (NUHMl) [58] the free parameters are 

Mo, Mh, Ml/2, ^0, tan/3, sgn(/i) . (1.2) 

In addition to the CMSSM parameters there is a universal scalar Higgs mass parameter 
at the unification scale, Mh^^ = Mh^ = Mh- It can differ from the other supersymmetric 
scalar masses at the unification scale. Mo. In NUHMl fits prior to the LHC data, light 
supersymmetric masses well below a TeV are also preferred [53-55, 59]. 

In 2011, the LHC experiments ATLAS and CMS have each taken approximately 5 fb~^ 
of proton-proton collision data at y/s = 7 TeV. The search for SUSY final states within these 
data has not yet revealed any signal. Instead, rather stringent bounds on the parameters 

^Employing the same minimal particle content but guaranteeing proton stability instead by baryon 
triality, there are more than 200 free parameters [46, 47]. 

■^Models with 19-21 new parameters were for example considered in [48-50]. 
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of the CMSSM are published in prehminary analyses [60]. The strictest bounds come from 
inclusive searches for events with missing transverse energy and jets in the final state. 

In this paper, we investigate in detail the consequences of these exclusions for the 
global interpretation of all existing data using our framework Fittino [61]. Fittino con- 
structs a global variable as a function of the model parameters and performs a Markov- 
Chain-Monte-Carlo (MCMC) scan of the parameter space. For the work presented here, 
138 million points in the parameter space have been found with a below 30. 

In analyzing these parameter scans we address the following questions: 

• What is the still allowed region in the CMSSM and NUHMl parameter space after 
two years of LHC data? 

• To what extent are the interpretation of the non-LHC measurements within the 
CMSSM/NUHMl and the non-observation of SUSY at the LHC in mutual tension? 

• What would be the impact of a light SUSY Higgs boson with a mass of around 
126 GeV, compatible with the latest LHC results [62, 63], on the CMSSM/NUHMl 
fits, and what would be the implication for its branching fractions and ratios of 
branching fractions? 

• What are the implications of the CMSSM/NUHMl best fits for direct and indirect 
searches for WIMP dark matter? 

• To what extent is the remaining preferred parameter region in the CMSSM/NUHMl 
fine-tuned? 

Technically, the methods employed to obtain the results presented in this paper differ 
in various aspects from our previous studies [56, 64, 65] and from those used in other 
global fits with LHC exclusions [66-76]. In contrast to our previous paper, the theoretical 
calculations no longer rely on the combination of codes in the Mastercode package [77] but 
make use of various individual codes for the different theoretical calculations (see Sec. 2). 
Particular emphasis is placed on an improved modeling of the LHC exclusions also in 
parameter regions away from the published 95% confidence level (CL) contours. This is 
achieved through a fast simulation of the signal for each tested point in parameter space 
instead of a simple parametrization of the CL. The experimental and theoretical input 
from direct and indirect WIMP searches is obtained from AstroFit [78]. In addition to 
the standard quantitative evaluation of fine-tuning [79] we develop a new, phenomenological 
method based on parameter correlations. 

The paper is organized as follows. In Sec. 2, the Fittino framework and the statistical 
methods to interpret the global variable are introduced. A comparison of frequentist 
and Bayesian interpretation of an example fit is also given. In Sec. 3, we present the 
experimental input employed in the fits. In particular, we consider indirect information 
from laboratory-based precision experiments, information from dark matter searches, and 
constraints from the exclusion of SUSY particles and Higgs bosons at previous colliders 
and at the LHC. We explain in some detail the method to estimate the LHC confidence 
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levels at arbitrary points in parameter space through fast simulation. The results of various 
fits are presented in Sec. 4. We first discuss global fits with updated observables excluding 
direct searches at the LHC, and then global fits including the current LHC exclusions. 
Then the impact of a potential Higgs signal is evaluated both within the CMSSM and the 
NUHMl models. This is followed by a discussion on the implications of the obtained fit 
results for direct and indirect dark matter searches. Subsequently, the fine-tuning of the 
best-fit regions is discussed. In addition, we explore parameter correlations as a measure 
of the required tuning of parameters with respect to each other. Furthermore, we present 
a new way to estimate theory systematics from a variation of the energy scale Q at which 
the weak scale mass parameters are evaluated after the RGE running. Our conclusions are 
given in Sec. 5. 

2 Obtaining constraints on model parameters 

In this section, we give an overview of the codes and methods used for the global fits of 
the CMSSM and NUHMl models. First, an overview over the codes used for scanning the 
parameter space and predicting the observables is given, followed by an explanation of the 
statistical methods used to derive the results. 

2.1 The Fittino framevi^ork 

The Fittino [61] framework is used to perform a global Markov Chain Monte Carlo 
(MCMC) scan (for earlier implementations and references see [56]) of the supersymmetric 
parameter space in all relevant parameter dimensions. For the calculation of the SUSY 
particle spectrum, the program SPheno [80] version 3.1.4 (SoftSUSY [81] version 3.1.7 as a 
cross-check) is interfaced. The resulting spectrum is used in micrOMEGAs [82] version 2.2 
for the prediction of the dark matter relic density, in FeynHiggs [83] version 2.8.2 for 
the prediction of the Higgs masses, the W boson mass my/, the effective weak mixing 
angle sin^ 9(,q, and the anomalous magnetic moment of the muon a^, in SuperlSO [84] ver- 
sion 3.1 for the flavor physics observables, and in AstroFit [78] for the evaluation of the 
direct and indirect detection of dark matter observables. A cross-check of the relic density 
computation between micrOMEGAs versions 2.2 and 2.4 [85] and DarkSUSY [86] version 5.0.5 
(via the AstroFit interface) is performed, which yields compatible results in terms of al- 
lowed parameter space. For the prediction of the branching fractions of the lightest MSSM 
Higgs boson, HDECAY [87] version 4.41 is used. The available limits on SM and non-SM 
Higgs bosons, including most available limits up to and including the ones presented by 
the LHC and Tevatron collaborations at the Spring Conferences 2011, are evaluated using 
HiggsBounds [88] version 3.2. The translation of the excluded CL into a contribution 
is explained in Section 3.3.2. 

The global is calculated in three steps: For all measurements Omeas given in Tab. 2, 
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- 4 - 



is calculated for each parameter point P where the sum runs over all A^i 
In case of upper bounds (e.g., the bound on 13{Bs — )• fj.fi)) 



meas 



measurements. 




2 



forO;,,,(P)>Of. 



limit 



(2.2) 



otherwise 



is calculated, where <t* is the assumed theoretical uncertainty of the prediction. Finally, 
the contributions from HiggsBounds, from the LHC SUSY search constraint, and from 
the direct and indirect detection of dark matter constraints are calculated, as outlined in 
Section 3, and added to the global x^- 

The MCMC algorithm employed uses a continuous optimisation of the width of the 
Gaussian proposal density functions (pdf) based on the variance of the accepted parameter 
points in all dimensions. Different settings for the ratio between the variance of recent 
accepted parameter points and the pdf widths are employed in parallel to combine fine 
scans of the minima with the ability to cover the complete parameter range. These methods 
and settings were developed specifically for the application in the fits presented here. They 
were found superior to our previous MCMC scan implementations [56, 64, 65], where the 
pdf was optimised in a separate study before performing the global MCMC fit. At least 
3 million points are obtained within Ax^ < 5.99 from the minimum for each individual fit. 

2.2 Statistical interpretation 

In the frequentist interpretation of the MCMC fit, first the point with the smallest x^ is 
identified. The 1-dimensional la (2-dimensional 2a) boundaries are defined by Ax^ < 1 
(Ax^ < 5.99) above the minimal x^ hidden dimensions are treated by the profiling 
technique, i.e. the hidden dimensions are scanned until the point with the lowest x^ ^or 
the given visible dimensions is found. The best fit point is the point with the smallest 
X^- Due to the excessive computing time necessary to find a reliable exact minimum, it 
is computationally prohibitive to perform toy fits [56], from which the exact CL coverage 
of the Ax^ < 1 and Ax^ < 5.99 contours could be derived reliably. Therefore, we cannot 
claim an exact match of the la uncertainty with 68 % CL. 

In the Bayesian interpretation, the full posterior pdf for all parameters is extracted from 
the MCMC local density. The m-dimensional marginalised posterior pdf for m parameters 
is drawn from the full n-dimensional posterior pdf by integrating out all other parameters. 



The corresponding la {2a) boundaries are defined by the smallest interval covering 68% 
(95%). The interval is built iteratively from the binned, marginalised posterior pdf by 
ordering all bins according to their probability pi, starting with the highest probability. 
Then the pi are subsequently added, until 68 % (95 %) is reached: 




(2.3) 
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Table 1. Results of the example fit for the frequentist and the Bayesian interpretation. For 
the frequentist interpretation the point with the smallest is given with the corresponding 1- 
dimensional la uncertainties. For the Bayesian interpretation, the global mode of the full posterior 
pdf in the case of flat priors is given, i.e. the point with the highest local 4D point density. In 
addition the maxima of the marginalised 1-dimensional posteriors arc shown with the boundaries 
of the smallest interval covering 68% around the maximum. 



parameter 


best fit (freq.) 


global mode (Bayesian) 


mode (marg., Bayesian) 


tan /3 

Ml/2 [GeV] 
Mo [GeV] 
^0 [GeV] 


10.3;9;5 

tro q+87.0 

<Jo-J-14.9 


8.2 
270.3 
52.4 
23.5 


q c:+21.3 

CO c+5055.5 
00.0_42 5 

4UO.U_2027.0 



The allowed region is the one which contains all bins in the range [1, maxnbin]- The remain- 
ing bins are outside of the allowed region. In this way, both the 2-dimensional 2a areas as 
well as the 1-dimensional (local) modes and the corresponding la intervals are constructed 
from the full posterior pdf. 

In Fig. 1, both interpretations are shown for an example fit in order to illustrate the 
differences. The fit is performed using the same setup as in [64]. Only the observables 
listed in Tab. 2 are used, in contrast to [64]. This example fit is only used for testing the 
statistical procedure, but not used for the results presented in Section 4. Figure 1(a) shows 
the 2-dimensional profile in the (Mo,Mx/2) plane. The small black contour on the left 
indicates the 2-dimensional 2a region in the frequentist interpretation. Figure 1(b) instead 
shows the MCMC scan point-density in the same plane. The black contour indicates the 
2-dimensional 95%-CL region for the same fit using a flat prior. In the white area, no 
points have been reached by the MCMC chain. Only scans with constant proposal density 
functions are used, in contrast to the continuous optimisation used for the other fits as 
described in Section 2.1 in order to keep the MCMC point density proportional to the 
posterior pdf. Note that for the frequentist interpretation also points which were rejected 
by the MCMC algorithm have been taken into account. This explains why there are 
small scanned areas in the {Mq, M1/2) plane which appear in the plot for the frequentist 
interpretation but not in the plot for the Bayesian interpretation. The results of the 
frequentist and Bayesian interpretation for the example fit are summarised in Tab. 1. 
Figures 1(c) and 1(d) show the Bayesian interpretation in the case of two different priors. 
Both priors have been chosen flat in Aq and tan/3, while the dependence on Mq and M1/2 
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Figure 1. Comparison of different statistical interpretations of the same example fit based on [64]. 
(a) frequcntist profile likelihood; (b) Bayesian marginalisation with flat prior; (c) and (d): two 
different non-flat priors, as described in the Text. The black lines indicate the 2a (frequentist) and 
95% CL (Bayesian) contours, respectively. 



is chosen to be 
P(Mo,Mi/2) 



oc 



Mo My 2 



[Fig. 1(c)], P(Mo,Mi/2) 



oc 



Mi + M2/^ 



[Fig. 1(d)] . 
(2.5) 



Given the significant prior dependence of the results, we employ the frequentist interpre- 
tation for the results given in this paper. In addition, the frequentist interpretation allows 
us to increase the point density of the MCMC scans by also including rejected points and 
by optimising the proposal density function continuously during the fit. The appearance 
of the so-called focus point region (allowed points at rather low M1/2 and very high Mq) 
in the Bayesian interpretation for the exactly same fit and the same MCMC chain as in 
the frequentist interpretation also helps to understand the significantly different allowed 
regions found in the literature in previous publications, see e.g. Fig. 1 in [57]. 
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Table 2. Low-energy observables employed. In general, experimental and theoretical uncertainties 
have been added in quadrature. The top quark mass mt is only used as an observable for the fit 
where mt is also floating as SM input parameter. 



B{b 57) 


(3.55 ±0.34) X 10"^ 


[89] 




< 4.5 X 10-9 


[90] 


B{B Tv) 


(1.67 ±0.39) X 10-"^ 


[89] 




17.78 ± 5.2 ps-i 


[89] 


exp _ SM 


(28.7 ±8.2) X 10-1° 


[5, 6] 


mw 


(80.385 ± 0.015) GeV 


[91] 


sin^ 9cs 


0.23113 ±0.00021 


[92] 




0.1123 ±0.0118 


[3] 


mt 


(173.2 ± 1.34) GeV 


[93] 



3 Experimental constraints 

In this section, we describe the present and potential future experimental measurements 
and searches which we employ to constrain the super symmetric parameter space. These 
can be grouped into three classes: 

1. Indirect constraints involving supersymmetric loop corrections; 

2. Constraints from astrophysical observations; 

3. Direct sparticle and Higgs boson search limits from high-energy colliders. 
We discuss these groups of observables in turn below. 

3.1 Indirect constraints 

Supersymmetric particles may contribute to various low-energy observables through loop 
corrections. The various types of indirect constraints that can be utilized to explore the 
SUSY parameter space are 

• rare decays of B-mesons; 

• the anomalous magnetic moment of the muon; 

• electroweak precision observables. 

The observables and their measured values as used in this paper are listed in Tab. 2. The 
different observables are briefly described in Sections 3.1.1 to 3.1.3. The Standard Model 
parameters that have been fixed are collected in Tab. 3. In one particular fit, the top quark 
mass is allowed to float. 
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Table 3. Standard Model parameters that have been fixed, apart from mt in one particular fit. 



1/aem 


128.962 


[6] 


Gf 


1.16637 X 10"^ 


[93] 


Q!s 


0.1176 


[93] 


mz 


91.1876 GeV 


[93] 


rub 


4.19 GeV 


[93] 


mt 


(173.2 ± 1.34) GeV 


[93] 


rrir 


1.77682 GeV 


[93] 


nic 


1.27 GeV 


[93] 



3.1.1 Rare decays of B mesons 

Strong constraints on new physics models can be obtained from flavor observables, in- 
cluding, in particular, B-meson decays. The precise measurement of branching fractions 
of rare decays which are helicity suppressed or which are mediated only at the loop level 
by heavy particles, places important restrictions on the supersymmetric parameter space. 
Significant constraints come from Bg oscillations, the branching fraction of B ^ tv and 
the inclusive branching fraction of 6 — )• 57. Recently there has been a substantial improve- 
ment in the limit on the branching ratio of the decay Bs — s- /i/i. The best limit is now 
13{Bs —7- fifi) < 4.5 X 10^^ [90], which is used as a default in our fits. In addition, we 
study the impact of a potential observation of Bg — t- fxfi by LHCb [94] at the SM rate of 
Ssm{Bs — )• fJ-n) = (3.2 lb 0.2) X 10~^ [95], i.e. assuming that no new physics is observed 
in this rare B-decay. Note that compared to our previous analyses [56, 64, 65], various 
observables from B- and K-meson decays have been discarded, as their contributions to 
the global fit have been found to be negligible. 

3.1.2 The anomalous magnetic moment of the muon 

Even though the anomalous magnetic moment of the electron is measured approximately 
two hundred times more precisely than that of the muon, (g — 2) = 2a^, the sensitivity 
of to new physics is enhanced by a factor of {m^/rrie)'^ ~ 43 000, and thus represents 
a much more powerful constraint on the SUSY parameter space. While the measurement 
of [5] is undisputed, there is an ongoing debate about the accuracy of its Standard 
Model prediction [6]. In particular, the non-perturbative contribution to from the 
hadronic vacuum polarization has to be extracted from experiment, via the data from 
e^e~ annihilation to hadrons or from r-lepton decays. For our base-line fit we have chosen 
the value of based on e+e~ data, which has been argued to be theoretically cleaner 
than that based on r decays. Note that in the CMSSM there is a strong correlation between 
B{Bs — )• /U/i) and [96]. 
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3.1.3 Electroweak observables 

In contrast to our previous analyses [56, 64, 65], the measurements of the Z-boson width 
and couphngs, and the hadronic cross-section on the Z pole have been excluded from the fit. 
Indeed, despite their high precision and the absence of any ambiguity in the interpretation 
of the measurement, the contributions of these observables to the fit have been found to 
be negligible. Relevant constraints only come from the effective weak mixing angle sin^ 9^^ 
and the W-boson mass. The values of the SM gauge couplings, and of the bottom, charm, 
and T-lepton masses are fixed inputs to the fit, see Tab. 3. 

3.2 Astrophysical constraints 

Further constraints arise from astrophysical and cosmological observations such as the cold 
dark matter relic density as well as from direct and indirect dark matter searches. We 
use AstroFit [78] for the evaluation of the direct and indirect detection of dark matter 
observables. The relic density is provided by micrOMEGAs [82]. AstroFit provides an 
interface to DarkSUSY [86] and adds an extensive database of relevant astrophysical data 
to the Fittino fit process. In particular, the particle spectrum and couplings of a given 
SUSY model are passed from Fittino using the SLHA standard [97, 98] to AstroFit which 
uses various DarkSUSY applications to calculate the predictions Othco for the requested 
observables O. This result is then compared to actual experimental data and contributes 
to the as Ax^ = [{Oexp — Otheo)/c'"exp]^; where Oexp is the experimental measurement 
and (Texp the la (combined statistical and systematic) error on it. For exclusion limits 
stated at a certain confidence level {e.g. from gamma-ray observations of dwarf spheroidal 
galaxies) or experimental results given in terms of confidence contours for claimed signals 
(as in the case of direct detection), AstroFit returns Ax^ = below the exclusion limit and 
inside the innermost confidence contour, respectively; outside this region, it interpolates 
between the A^^ values that correspond to the stated confidence levels and extrapolates 
quadratically beyond the largest explicitly given CLs (for limits, this is done by assuming 
that O follows a Gaussian distribution around zero). On return, AstroFit passes the Ax'^ 
contribution from every included observable back to Fittino. In this paper, we use the 
dark matter relic density as well as the neutralino scattering and annihilation rates that 
are relevant for direct and indirect matter searches. 

3.2.1 The cold dark matter relic density 

Evidence for a considerable cold dark matter component in the composition of the Uni- 
verse derives from a great number of observations that cover a large range of distance 
scales. In particular, the WMAP satellite has determined its cosmological abundance from 
observations of the cosmic microwave background to an impressive accuracy of J^dm^^ = 
0.1123ib0.0035 [3]. Weakly interacting massive particles (WIMPs), if neutral and stable on 
cosmological time scales, are very good candidates for the cosmological dark matter - the 
prime example being the lightest neutralino [99]. Thermally produced in the early Universe, 
the WIMP relic density is set by the (co-)annihilation rate {av)DM, see e.g. Refs. [100, 101]; 
requiring that all dark matter is constituted by thermally produced neutralinos makes the 
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relic density one of the most constraining observables for the CMSSM parameter space. 
In our fits, we calculate the relic density for comparison with both micrOMEGAs [82] and 
DarkSUSY [86] (via the AstroFit interface). We add a 10% theory uncertainty in the cal- 
culation of the relic density due to missing higher order corrections in the anihilation cross 
sections [102-104]. 

3.2.2 Direct dark matter detection 

Direct detection instruments look for elastic scattering events of WIMPs with nuclei of the 
according target material in underground laboratories. They detect signals via scintillation, 
phonons or ionization, depending on the experiment (see e.g. [105] for a review of direct 
dark matter detection). The spin-independent elastic scattering cross-section per nucleon, 
cJsi, conventionally adopted for a model- independent comparison between experimental re- 
sults, is calculated with AstroFit, taking properly into account the potentially different 
scattering cross-sections on protons and neutrons. Theoretical uncertainties of approx- 
imately 20% from the so-called sigma-term [106] affect the calculation through changes 
to the form factor [105]. These uncertainties are estimated by varying the form factor 
within DarkSUSY and evaluating the resulting spread of the predicted value of csi for dif- 
ferent points in the CMSSM parameter space along the expected exclusion. A resulting 
systematic uncertainty on cJsi of up to 50 % was found and added in quadrature to the 
experimental uncertainty obtained from the exclusion contour as described above. 

Recently, controversial results have been reported from different direct detection ex- 
periments. While some of the collaborations, DAMA/LIBRA [107], CoGeNT [41, 108] and 
CRESST [109], promulgate the detection of a signal, others, primarily XENONIOO [39], 
have shown upper limits nominally excluding the regions of signal detection. We adopt 
the current XENON limit as well as various projections including the XENONIT exper- 
iment as given in Ref. [110]. For an upper limit such as from XENON, the application 
of the aforementioned theoretical uncertainty on cjsi of 50 % translates into a maximally 
achievable Ax^ contribution of 4 for the XENON limit. For the positive signals claimed 
by DAMA/LIBRA, CoGeNT and CRESST, no such hmit applies directly. For those, we 
observe prohibitively high Ax^ contributions and thus do not provide further fit results in 
Section 4. 

3.2.3 Indirect dark matter detection 

Indirect detection instruments - tracing products from dark matter annihilation such as 
photons, antiprotons, positrons, and neutrinos in the corresponding cosmic-ray fluxes - 
have been included via the gamma-ray channel (see e.g. [Ill] for a review of indirect de- 
tection methods). Upper gamma-ray flux limits from dwarf spheroidal galaxies, preferred 
observation targets due to their high mass-to-light ratio and low background, from obser- 
vations with Fermi-LAT [112] and H.E.S.S. [113] have been included in the fit through 
AstroFit [78]. 

We note that gamma-ray limits from dwarf spheroidals lead to the currently strongest 
model-independent bounds on the annihilation rate (at least for the range of dark matter 
masses that is relevant here). In this work, we therefore do not include available photon 
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flux information from other sky regions, such as the galactic halo. The same holds true for 
other indirect detection channels, in particular antiprotons and positrons. 

3.3 Direct search limits from high-energy colliders 

Searches for supersymmetric particles have been performed at various high-energy colliders, 
with strong constraints on the MSSM particle spectrum, including supersymmetric Higgs 
bosons, specifically from LEP, the Tevatron, and the LHC. In the fits we include the limit 
on the lightest chargino mass from LEP, limits on the supersymmetric neutral and charged 
Higgs boson masses from LEP and the Tevatron, and, in particular, limits from recent 
LHC searches for supersymmetry in channels with jets and missing transverse energy. 

3.3.1 The limit on the lightest chargino mass from LEP 

The LEP collaborations have searched for charginos within the MSSM. The analyses assume 
gaugino mass unification at the GUT scale. Mi = M2 = M3 = which leads to the 

relation Mi = | tan^ 0^yM2 at the electroweak scale. This holds for both the CMSSM 
and the NUHMl considered here. The LEP collaborations searched for the dominant 
decay mode Xi ~^ l^^*Xi [H^]. We have used the most conservative chargino mass 
limit of m ± > 102.5 GeV, including a conservative estimate of the theoretical uncertainty 

of Am ± = 1 GeV. Because of gaugino mass unification, the limit on the chargino mass 

Xi 

indirectly also excludes neutralinos with masses m^o < 50 GeV in the constrained models 
considered here [93, 115]. 

3.3.2 Limits on Higgs boson masses 

The program HiggsBounds [88] is used to implement the limits on the supersymmetric neu- 
tral and charged Higgs boson masses from LEP, the Tevatron and the LHG. HiggsBounds 
first determines the strongest search channel for the various Higgs signatures at every tested 
model point, based on the expected limit. For the SUSY models investigated here, the LEP 
limits on neutral Higgs bosons prevail in most of the relevant parameter space. The CLs+t 
and Chg values corresponding to these limits are tabulated in HiggsBounds for all signal 
strengths. Thus, even though the likelihood of the data for a given Higgs model is not 
directly available, the Xh contribution from the Higgs searches can be calculated from the 
observed CLg+b value for the most sensitive search channel and for any given Higgs mass 
and signal strength, with the valid assumption of a Gaussian distributions of CLg+b- The 
relation 

X^ = 2[erf-i(l-2CL,+fe)]2 (3.1) 

is used to determine the full x| contribution disregarding theoretical uncertainties of rrih. 
The contribution of the Higgs limits to the total of the fit is then obtained from the 
x|(m./j) distribution after folding with a Gaussian of width Am/j = 3 GeV to take into 
account the theory error [116-118] of the Higgs mass calculation. This procedure is imple- 
mented as an addition to HiggsBounds [119, 120]. 

In addition to the LEP results discussed above, the LHC collaborations have now 
presented an updated result on the combination of searches for the Standard Model Higgs 
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boson [62, 63] with a combined allowed range at the 95% CL limit of 117.5 < mh < 
127.5 GeV, exceeding the LEP limit. No detailed translation of the CLs+b values from the 
LHC-experiments into contributions is possible given the published values. Since the 
lightest Higgs boson in the CMSSM is SM-like, we can nevertheless study the effect of the 
new limit on the global fit in an approximate way by cutting on irih, taking into account 
also the theoretical uncertainty Am/i = 3 GeV. Furthermore, in a separate study we discuss 
the impact of a hypothetical SM-like Higgs boson signal at around ruh = 126 GeV [62, 63] 
on global fits in the MSSM and the NUHMl. 

3.3.3 SUSY searches at the LHC 

At the LHC, the most stringent limits on supersymmetric models with R-parity conser- 
vation are obtained from searches in channels with jets and missing transverse energy, 
EJp^^^ [60] , where squarks q and gluinos g are produced in pairs and decay through q — t- qxi 
and g — t- qqxi to purely hadronic final states with missing energy from the weakly inter- 
acting and stable lightest neutralino Xi- 

In order to properly include the SUSY exclusions from the LHC searches in the global 
SUSY fit, it is not sufficient to only consider the 95% CL bounds published by the exper- 
imental collaborations for specific models and particular choices of parameters. Instead, 
it is mandatory to determine the corresponding x^ contribution of the LHC observables, 
Xlhc every point in parameter space. To that end we emulate the experimental SUSY 
search analyses by using state-of-the-art simulation tools and a public detector simulation 
to approximate the signal expectation for the specific model under consideration and any 
given set of model parameters. Details of the analysis are described below. 

We determine the LHC SUSY exclusions from the search for squarks and gluinos using 
final states with 2, 3 or 4 jets, zero leptons, and missing transverse momentum, following 
a recent ATLAS analysis [121]. The signal cross-section is simulated on a grid in the 
{Mq, Mi^2) plane, fixing tan/3 = 10, Aq = and /i > 0. (We discuss the dependence 
of the signal on tan/3 and Aq below.) First, for each point of the (/Vfoi-^1/2) grid, the 
SUSY mass spectrum is calculated with SPheno 3. 1.0 [80]. Events are then generated 
using Herwig 2.4.2 [122] and passed on to the fast detector simulation DELPHES .1.9 [123] 
where reconstruction and detector smearing are performed. The efficiencies and fake rates 
of particle identification given by DELPHES are corrected using public measurements of 
detector resolutions by ATLAS [124]. Finally, the events have to pass the selection cuts 
specified in [121]. 

The signal estimate is normalized to the NLO-I-NLL SUSY-QCD prediction for the 
inclusive squark and gluino cross-sections [125-131] and the NLO prediction for the pro- 
duction of electroweak sparticles [132]. While the theoretical uncertainty of the signal 
cross-section varies along the parameter space, we assign a constant 30 % error to the sig- 
nal yield, which is a conservative estimate in the vicinity of the current LHC sparticle 
exclusions. The variation of the signal yield with the remaining CMSSM parameters tan /? 
and Aq has been studied carefully. We found that it is always well within the systematic 
uncertainty, as exemplified in Fig. 2 for two points in the (Mq, Mx/2) parameter space. 
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Figure 2. Simulated signal yield for events with 2, 3 or 4 jets for two points in the {Mq, M1/2) 
parameter space. Shown are the SM background (gray), the CMSSM signal expectation for tan/3 ~ 
10 and Aq = (yellow) together with the systematic uncertainty of 30% (orange). Also shown 
through different markers are various other central value signal estimates based on different values 
of tan P and Aq . 



Our estimate of the SM background contribution is based on the ATLAS analysis [121], 
which reports the results and statistical interpretation of the ATLAS fully hadronic inclu- 
sive SUSY searches using 165 pb~^ of integrated luminosity. A systematic uncertainty 
of 23% on the background has been estimated in [121]. Using the ATLAS background 
and background systematic uncertainty, together with a systematic uncertainty of 30 % 
on the signal cross-section, we can exactly reproduce the published ATLAS limits at 
£'■1* ^ i65pb'i and = Ifb-^ The ATLAS analysis for = 4.7fb-\ however, 
uses a profile likelihood based interpretation and six signal regions with tighter cuts. Be- 
cause of the large amount of background control regions, this method cannot easily be 
implemented in an exact way within the fast Monte Carlo analysis used here. Most impor- 
tantly, the background systematics are controlled within the profile likelihood technique. 
Therefore, we approximate the generic LHC results by reducing the systematic uncertain- 
ties used in the likelihood ratio analysis within realistic bounds and by slightly increasing 
the signal efficiency. Fig. 3 shows that this fast Monte Carlo technique is able to describe 
the LHC results within the 1 a variation of the expected limit, both in scale as well as in 
shape. Thus, this method should provide a good description of the Xlhc contribution at 
every point in SUSY parameter space. 

The exclusion limit is calculated in the parameter space for each point with expected 
yields for signal s and background b and observed yield n, calculating a test statistic 
t = —2lnQ, with Q being the likelihood ratio: 

L(^_ (3.2) 

L[b, n) 

with the Poisson statistics L{ii,n) = /i'^exp"^ /n\, the number of expected and observed 
events, fj, and n, respectively, for the given point. The value of s is a function of the SUSY 
parameters, whilst b is fixed. The additional uncertainties are taken into account by a 
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Figure 3. contribution from the LHC SUSY search implementation compared to the pubUshcd 
ATLAS and CMS Hmits [60] . Good agreement of the estimated hmit with the expected hmit of the 
LHC coUaborations is achieved. 

smearing of the Poisson distribution, and we consider a signal excluded with 95 % CL if 

CL,+f,= / P,+fe(t)di <0.05, (3.3) 

^obs 

with -Ps+fe the probability density function of t assuming the presence of a signal, and ^obs 
the actually observed value of t. A given CLs_|_ft value can be approximately translated into 
a contribution using [120]: 

x' = 2[erf-i(l-2CL,+,)]2. (3.4) 



To obtain the expected exclusion limits the Asimov data set n = b was used. 



4 Results 

In this section, we present the results of our CMSSM and NUHMl fits for various sets of 
input observables. We first summarize our general findings and then discuss the features 
of the different fits in more detail in Sections 4.1 to 4.5. 

In Tables 4 and 5, we display selected results of the fits within the CMSSM and the 
NUHMl for various sets of input observables. For all fits we require the lightest neutralino 
to be the LSP, consistent radiative electroweak symmetry breaking and the absence of 
tachyons. 

In total, 7 different input parameter sets have been tested for the CMSSM. The fits 
based on the CoGeNT^ [108] direct dark matter detection claim yield a prohibitively 

^The fit involving DAMA/LIBRA data has been omitted entirely since the CoGeNT regions lie closer 
to the accessible space in the CMSSM. 
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Table 4. Summary of the results for various CMSSM fits with different sets of input observables. 
The names of the fits and the observables included are explained in the text. 



Fit 


Mo [GeV] 


Ml/2 [GeV] 


tan P 


^0 


X^/ndf 


LEO 
LHC 
LHC+XENONIT 
LRC+nih = 126 GeV 


296.ltJ366.8 


375.41^^^/ 
74 7 4+303.4 

1167 4+594.0 


14.91^^2^ 

34.4;^?i 

98 q+21.2 

oq 0+16.7 


186.3li^:^ 

o 7+5266.3 
OiO. /„2166.1 

-2969.ll?il:8 


10.3/8 
13.1/9 
15.0/9 
18.4/9 



Table 5. Summary of the results for the NUHMl fit. 

Fit Mo [GeV] Ml/2 [GeV] M|f [GeV^] tan/3 Aq X^/ndf 

NUHMl+TO,, UA.itlli 655.51^^*0° (-l-^l":?) x 10^ 29.4l'?j -511.21^^^;^ 15.3/8 

large minimal x^, implying that these potential DM observations cannot be accommo- 
dated within the CMSSM. We thus do not discuss them further. For all fits we require 
the lightest neutralino to be the LSP, consistent radiative electroweak symmetry breaking 
and the absence of tachyons. Note that a positive sign of /i is preferred to describe the 
anomalous magnetic moment of the muon, so we have fixed sgn(/_f) = + in our fits. 

The fit labeled "LEO" ("low energy observables") in Tab. 4 is performed excluding 
the direct search for SUSY at the LHC, but including all constraints from supersymmetric 
loop corrections (Section 3.1), astrophysical observations (3.2) and sparticle and Higgs 
boson search limits from LEP and the Tevatron (3.3.1 and 3.3.2). The best fit point 
and the 1-dimensional la parameter uncertainties are in excellent agreement with earlier 
studies [64], despite the fact that only a fraction of the observable set from [64] is used 
here. As mentioned in Sections 3.1.1 and 3.1.3 various flavor and electroweak precision 
observables have been omitted from the fit, as they do not add significant information for 
constrained SUSY models. Instead they add a constant offset to the x^j which is nearly 
independent of the SUSY parameter point or model under study. The offset only depends 
on the SM parameters and the SM precision observables themselves. In Ref. [56] it has 
been shown that the SM parameters decouple completely from the CMSSM parameters for 
the observables used here. Therefore, a reduced observable set results in a clearer message 
on the actual agreement of the SUSY prediction with the relevant data for a given model 
point. This is shown in Fig. 4. The red lines show the profile of the x^ minimum of the fit 
as a function of Mq, once with the full observable set from [64] (dark red), and once with 
only the set of observables described in Section 3 (light red), but using the theory codes 
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Figure 4. Comparison of fits witli small/large set of observables and witli and witliout exact 
HiggsBounds Xh- The minimum of all fits apart from the one using HiggsBounds for the calculation 
of the full Higgs contribution is set to Ax^ = for better comparison. Wc show here only the 
dependence on Mq. M1/2, and tan/3 are profiled out. 

and observable vales from [64]. No significant difference can be observed. 

An interesting effect can be observed for the observable B{B — )• tv): While SUSY 
points with strong disagreement to the measurement exist, the majority of the allowed 
region exhibits a prediction of B{B — )• ti') near its SM value, since the SUSY contribution 
tends to decrease the branching fraction, leading to stronger disagreement than in the SM. 
Thus, also B{B — )• tu) contributes to a constant offset of the with respect to previous 
results with lower observed values of the branching fraction (e.g. [56]), while having no 
significant effect on the allowed region above the minimum. 

In determining the number of degrees of freedom, ndf, we include one-sided bounds as 
measurements, also if the minimum of the fit fulfills the bound. This is justified because 
the model exhibits parameter points violating the bounds. Removing the bound from the 
ndf count as soon as the minimum or the regions directly around it satisfy it, would mean 
that the model is penalized for being consistent with the data on the ground of constraints 
from other observables. Therefore, all bounds with potential contribution to the ^■re 
counted in ndf. 

Using the new observable set, and the new set of codes for the prediction of the 
observables as outlined in Sections 2 and 3, we obtain a x^ /ndf of 10.3/8, hinting at a worse 
agreement than using the full observable set [64]. However, as outlined in Section 2.2, no 
exact P-value can be calculated due to computational limitations and the non-Gaussian 
nature of the fit. Clearly the overall fit quality of the CMSSM is reduced when using 
the reduced observable set of Tab. 2, and taking the improvements in the sensitivity of 
the measurements between the old and new observable set into account. Additional small 
effects stem from the update to newer versions of theory codes as outlined in Section 2.1. 
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Another significant difference between the fit results from [64] and the results presented 
here is the more precise treatment of the LEP Higgs limits. In [64], HiggsBounds was used 
to calculate the exact position of the 95 % CL on nifi for each model point. Then, the left 
half of a parabola corresponding to a theoretical uncertainty of Arrih = 3 GeV and with 
the minimum oi Xh = ^ at the exact 95 % CL point was used to calculate Xh- Model points 
above the actual limit received a contribution of Xh — 0) while model points below received 
a Xh contribution corresponding to the theoretical uncertainty. 

For the results presented here, HiggsBounds is used to calculate the exact Xh 
every rrifi below and above the 95 % CL bound, folded with the theoretical uncertainty 
(see Section 3.3.2). Since for the best fit point of the LEO input set we obtain m/, ~ 
(113 lb 3) GeV, adding the actual Xh contribution to the minimum makes a significant 
difference in the overall x^ of the minimum. This can be seen from the blue lines of Fig. 4. 
The minimum in the dark blue curve is lifted significantly with respect to the fits without 
HiggsBounds (light blue). Since larger values of Mq > 2TeV yield ruh > 118 GeV, Xh 
does not contribute to the fit for large Mq and the two lines are in perfect agreement for 
large Mq. Thus, a thin stripe of the focus point region at large Mq and very small M1/2 
is included in the allowed region, see Sec. 4.1. Its width in M1/2 is mostly constrained 
by requiring m + > 102.5 GeV. The pull plot in Fig. 5(a) shows the contribution of the 
individual observables to the best fit point of the LEO fit. 

The fit labeled "LHC" in Tab. 4 uses the same input as the LEO fit, plus the results 
from the LHC SUSY search described in Section 3.3.3. The effect of the LHC limit is a shift 
of the minimal x^ and an increase of the number of degrees of freedom to Xi^ /ndf = 13.1/9. 
This clearly shows a further reduced agreement of the fit with the data, and the minimum 
is found at higher Mq, and tan/3. 

From the pull plot in Fig. 5(b) one can deduce that the LHC limit pushes the fit 
into a region where, in particular, is not described very well anymore. In comparison, 
the direct contribution of Xlhc small at the best fit point. The reason is that the x^ 
contribution from the low energy observables is almost constant for large enough Mq and 
Ml/2- Thus, if the LHC pushes the fit above a certain threshold in the mass scale, the 
minimum can also move higher until Xlhc ^'-'^ dominant anymore. In addition, x^ also 
prefers higher sparticle masses, balancing the push to lower mass scales from the precision 
observables. A further consequence of the LHC exclusions is that the preferred Higgs mass 
is shifted upwards to ruh 117 GeV. 

The fit labeled "LHC+XENONIT" in Tab. 4 uses all the results from Tab. 2 and the 
direct LHC SUSY search, but the current XENONIOO limit replaced by the projection 
for XENON IT as outlined in Section 3.2.2, assuming no observation. The XENON IT 
projection increases the minimal x^ to 15.0. In order to achieve lower direct detection 
cross-sections, the fit is pushed to somewhat lower values of tan /3 and excludes large 
positive values of Aq, resulting in an Aq dependence which is almost symmetric around 
^0 = 0. The detailed reasons are given in Section 4.4. 

We have also considered the possibility that B{Bs — )■ /i/u) is measured at its SM value 
oiB{Bs Hfi) = (3.2 ±0.3) x 10"^ at LHCb [94] (instead of the otherwise employed limit 
of 13{Bs — T- /i/u) < 4.5 X 10~^). Since the minimum of the LHC fit predicts B{Bs — t- /x/i) 
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Figure 5. Contribution of the individual observables listed in Tab. 2 to selected fits, (a) the LEO 
fit; (b) the LHC fit; (c) the LHC+m,i=f26 CMSSM fit; and (d) the LHC+m,, =126 NUHMl fit. 
In addition the limits from indirect WIMP detection and the chargino limit at LEP count for the 
calculation of the number of d.o.f. 



near its SM value, no significant change in fit probability or best fit point is observed, see 
the discussion in Section 4.2. 

The new limits on the SM Higgs boson mass of 117.5 < ruh < 127.5 GeV at 95 % CL [62, 
63] from the ATLAS and CMS experiments are not included in the detailed calculation in 
HiggsBounds. Therefore, no detailed calculation of the contribution from that constraint 
is performed. However, applying a simple cut of 114.5 < nifi < 130.5 GeV (including a 
theoretical uncertainty of Am/j = ±3 GeV) on the LHC fit does not lead to a significant 
effect on the allowed parameter space, since for the largest part of the parameter space in 
the LHC fit ruh is predicted in the allowed range due to the heavy sparticle mass scales. 

Finally, we explore a possible Higgs observation at rrih = (126ib2(exp) ib3(theo)) GeV, 
as outlined in Section 3.3.2. The corresponding fit is called "LHC+m^ =126" in Table 4. 
Apart from exchanging the HiggsBounds result for Xh for the assumption of a direct mea- 
surement, the inputs to the fit are identical to the LHC fit. The high sparticle masses 
required to yield such a high Higgs mass push the best fit point into the focus point region 
at Mo > 4TeV. The x^/ndf = 18.4/9 shows that ruh > 125 GeV is hardly compatible 
with a highly constrained SUSY model such as the CMSSM. This is also obvious from the 
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Figure 6. Parameter distributions for the LEO fit. Details about the best-fit point, marked by a 
star, are given in Table 4. In (a), the profile likelihood projection of the Ax^ < I (I-dimensional 
Icr) [red] and Ax^ < 5.99 (2-dimensional la) [blue] regions of the full 4-dimensional fit around the 
minimum are shown into the (Mo,M2/2) plane. In (b), the same is shown in the (ylo,tan/3) plane. 



pull contributions in Fig. 5(c). While a sufficiently high value of rrih = 124.4 GeV can be 
reached for the best fit point, the necessary high sparticle mass scales push the fit into a 
region which is incompatible with (^f — 2)^ and S-physics observables at the 2 to 3 o" level. 

By decoupling the Higgs sector from the squark and slepton sector in the NUHMl 
model, this problem can be partly remedied. As shown in Tab. 5 and Fig. 5(d), a slightly 
better fit is possible for the NUHMl {x^/ndf = 15.3/8) than for the CMSSM fit with 
m/i = (126 ± 2 lb 3) GeV, albeit with reduced ndf. Also, the added degree of freedom 
deepens the profile for Mq and around the minimum, since the Higgs limit does 

not push the fit into an area with limited agreement between predictions and indirect 
observables. For a more detailed description of the NUHMl fit, see Section 4.3. 

In the following sections, the results for individual input observable sets are discussed 
in detail. 

4.1 Fit without LHC exclusions 

In Fig. 6 the resulting parameter regions are shown for the LEO fit. The best-fit point is 
marked by a star. All hidden dimensions are profiled. As in all further plots of the allowed 
parameter space, the allowed regions display the full Ax^ < 1 (1-dimensional lo") range 
in red and the Ax^ < 5.99 (2-dimensional 2a) range in blue. The 1-dimensional la range 
is shown in order to allow reading off the single parameter uncertainties directly from the 
projection of the red areas onto the horizontal and vertical axes. 

The resulting plot does not represent a parameter scan in the shown dimensions with 
fixed hidden dimensions, but the full uncertainty regions for the global fit of the four 
CMSSM parameters. As shown in [64], including the SM parameters in the fit does not 
change the result as long as the presently available observables are used. This may change 
in the future, e.g. once very precise measurements of nih or sparticle masses are available. 
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(a) (b) 

Figure 7. Best-fit regions in the CMSSM including LHC results. Details about the LHC best-fit 
points are given in Table 4. Results for the la and 2(T regions are shown in comparison between 
the LEO and LHC fits in (a) the (Mq, M1/2) plane and in (b) the {Aq, tan /3) plane. The parameter 
projections and uncertainties are as in Fig. 6. 



In contrast to the results from [64] , and in contrast to the profile likelihood result from 
Fig. 1 which uses the observable set and prediction codes from [64], the focus point region 
at low Ml/2 ™d high Mq is allowed within the 2-dimensional 2a uncertainty range. This 
is explained by the cut through the profile already shown in Fig. 4. The focus point 
region is constrained to M1/2 < 200 GeV by the presence of the cut on m^± > 102.5 GeV 
described in Section 3.3.1. The minimal shape in Fig. 4 also explains why the 2a allowed 
region is so much larger than what would be expected from the la region. The fit exhibits a 
rather narrow minimum, where a good agreement with {g — 2)^ and JIdm can be achieved. 
For larger values of Mq and Mx/2) there is an almost constant disagreement with both 
measurements at the 2 — 3(T level. Hence, no further discrimination is achieved and the x? 
profile becomes almost flat, still within the < 5.99 range above the minimum. 

4.2 Fits with LHC exclusions 

In this section, we discuss the allowed CMSSM parameter space for different inputs from 
the LHC. The inputs with the strongest impact are the inclusive direct searches for SUSY 
at the ATLAS experiment described in detail in Section 3.3.3 and the search for Bg — )• Hjji 
described in Section 3.1.1. The additional very strong constraint stemming from a possible 
measurement of the lightest Higgs boson mass is studied separately in Section 4.3. Fig. 7 
shows the allowed parameter range of the LHC fit introduced in Table 4. The difference 
between the LHC and the LEO fit is significant. The position of the best-fit point is lifted 
outside of the directly accessible range of sparticle searches at low Mq and Mi 12 ~ 650 GeV. 
The focus point region is excluded by the LHC SUSY search up to our maximum value 
of Mq = 3.5 TeV of the LHC limit implementation. Also the LHC collaborations do not 
publish search results for larger values of Mq, but it can be assumed that a dedicated 
interpretation at LHC would exclude large parts of the rest of the focus point region from 
the fit (i.e. the two islands above Mq > 3.5 TeV in Fig. 7(a)). Due to the diminishing 
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(a) (b) 

Figure 8. Parameter distributions for the LHC fit {B{Bs — > < 4.5 x fO^^) and the 

LHC^''"(^=~'^^' fit in the (Mq, A/1/2) in (a) and in the (Ao,tan/3) plane in (b). Details about 
the best-fit points are given in Table 4. 



relevance of large SUSY mass scales for the solution of the hierarchy problem and EWSB, 
this is not further pursued here. 

It can be observed in Fig. 7(b) that the best fit point of tan f3 is shifted upwards 
significantly due to the direct LHC SUSY search limits, which do not affect tan (3 directly, 
as shown in detail in Section 3.3.3. This is the result of the correlation between tan/3 on 
the one hand and M1/2 on the other hand due to JIdm and {g — 2)^. A decent agreement 
in those variables for larger M1/2 implies higher scales of tan/3. 

In Fig. 8 we use the LHC limits as outlined above and in addition study variations 
of the limit on 13{Bs — )■ ^/u). In Fig. 8(a) and 8(b), the results of the fit including a 
potential measurement at the SM value I3{Bs — s- fifi) = (3.2 it 0.3) x 10^^ (assuming 
a small theoretical and experimental error of ±0.3 x 10~^ for the sake of studying the 
strongest possible impact) at LHCb [94] and possibly ATLAS and CMS is shown. Since 
the best-fit point of the LHC fit predicts B{Bs — )• ///u) close to its SM value, this has only a 
marginal effect on the minimal of the fit and the allowed parameter space. Note that the 
situation would have been somewhat different for a potential observation of B{Bs — )• 1^1 fi) 
at CDF [133], which would have clearly disfavored the focus point region. 

We show the predicted range of sparticle and Higgs boson masses for the LHC fit 
in Fig. 9(b) compared to the same prediction for the LEO fit in Fig. 9(a). As discussed 
above and in Fig. 4, the flat profile for large Mq and M1/2 predicts sparticle masses 
at m ~ 1 TeV for both the LEO and LHC fits. The only particle with a very strong 
upper bound on its mass is the lightest Higgs boson. The prediction of 711^, changes from 
rrih = (113.5 ± 3) GeV for the LEO fit to ruh = (117 ± 1.5) GeV for the LHC fit. The 
1-dimensional la range quoted here is clearly incompatible with the potential signal at 
TUfi = 126 GeV. However the profile is flatter for higher rrih, making it marginally 
possible to reconcile a heavier Higgs mass with the CMSSM, as outlined in Section 4.3. In 
addition, the lightest neutralino is bound from above at m^o < 500 GeV, as a result of the 
combined constraints on J^dm and (g — 2)^. 
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(a) (b) 

Figure 9. Predicted distribution of spartiele and Higgs boson masses from (a) the LEO fit and (b) 
the LHC fit. The full uncertainty band gives the 1-dimensional 2(t uncertainty of each mass defined 
by the region Ax^ < 4 after profiling over all hidden dimensions. 

We note that also in the LHC fit there is still enough room for slepton and gaugino 
masses below m < 500 GeV, which would be observable at an e^e~ linear collider with 
\fs < ITeV, despite the stringent search limits from the LHC. This is partly due to the 
fact that the LHC allowed 2a region extends slightly below the 95 % CL exclusion. The 
increasing LHC contribution for smaller spartiele mass scales is partly canceled by the 
decreasing contribution from the LEO observables. The stronger effect however is the fact 
that even highly constrained models as the CMSSM still exhibit enough freedom in the 
parameter space to reconcile light gauginos and sleptons with heavy squarks and gluinos, 
as shown in Fig. 9(b). This feature would be even more pronounced in more general SUSY 
models, in which the colored and the non-colored sectors explicitly decouple. 

4.3 Adding a potential Higgs boson mass measurement 

A very strong constraint on the CMSSM is expected from a direct measurement of the 
lightest Higgs boson mass rah,. This is shown explicitly in Fig. 10, where we display the 
profile of m/j for different fits: the LHC fit, the LHC+m/j =126 fit with fixed m^, the same 
with rrit = 173.2 it 1.34 GeV floating, and the NUHMl fit with the same observables. It can 
be seen that the CMSSM fits (red and green curves) can barely incorporate ruh = 126 GeV. 
At this edge the x^ rises dramatically excluding any higher Higgs boson mass values. This 
limit is just slightly increased by fioating mt (green curve), but the overall /ndf of this 
fit is still at the same rather inconsistent level of 18/9 as for the fit with fixed mt, for which 
18.4/9 is observed. It is obvious that the NUHMl fit, which decouples the Higgs boson 
mass parameter squared = Mfj^ = Mfj^ from Mq at the GUT scale, shows both a 
significantly better agreement of /ndf = 15.3/8 and the ability to accommodate higher 
values of m/,. 

In Ref. [96] a strong correlation between B{Bs — >• /i^) and was observed for fixed rrih. 
It was shown that the correlated region for the current experimental values of B{Bs — )• /i^) 
and is inconsistent with a Higgs boson mass at ruh = (126 ± 2(exp) it 3(theo)) GeV. 
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Figure 10. The dependence of the minimal of the fit on nih for different input observable sets and 
for the CMSSM and NUHMl. CMSSM and NUHMl fits with and without mn = (f 26 ± 2 ± 3) GeV 
are shown. For the CMSSM, whose best-fit-point coincides with a prediction of the SM value of 
B{Bs — >■ ^/x), the strongest impact comes from including nih = (126 ± 2 ± 3) GeV into the fit 
(compare dark and light red lines). Adding nit as a fit parameter has limited effect (compare green 
and light red line). For the NUHMl, a heavier Higgs boson is in principle easy to accommodate, 
however a heavy Higgs boson is disfavored in the NUHMl by the present bound on B{Bs fi^J.) [90] 
compared to a previous looser limit [134] (compare light and dark blue lines). 

Accordingly the pull plot Fig. 5(c) has a poor fit for a^. Performing a fit just for the 
three observables B{Bs — )■ /i/x), a^, and ruh = (126 ± 2(exp) it 3(theo)) GeV results in 
a Ax^ penalty of 7 with respect to just fitting B{Bs — )■ /X/u) and a^. This explains the 
strong difference between the fits with and without rrih = (126 it 2(exp) it 3(theo)) GeV. 
To further exhibit the combined effect of the limit on B{Bs — t- fifj.) and a potential Higgs 
boson observation, we show in Fig. 10 the profile of against m/j for different options 
of the measurement of B{Bs — t- fifi). While a measurement at the level of the SM of 
B{Bs —7- fifj.) = (3.2 lb 0.3) X 10~^ has no strong effect on the maximal achievable value 
of m/j for the CMSSM fit, there is a significant difference for the NUHMl fit, where large 
nifi is disfavored for a low B{Bs — s- ^/u). Even the current LHCb limit severely impacts 
the fit quality with respect to a looser bound on B{Bs — s- fifJ.), since the NUHMl fit with 
nih = (126±2(exp) ib3(theo)) GeV prefers B{Bs — fijj-) > 6.5 x 10~^. Of course this is not 
a general observation, since the i?-physics and the Higgs sector observables are not directly 
governed by the same parameters in the MSSM, while they are connected via GUT scale 
assumptions in the CMSSM and NUHMl. Therefore, also in this example, testing more 
general SUSY models than the CMSSM and the NUHMl would be beneficial. However, 
the current limitations on the correct implementation of the LHC results, as outlined in 
Section 3.3.3, has to be overcome. 

Figure 11 shows the allowed parameter range in the (Mo,Mi2) and (tan/3, Aq) planes 
for the CMSSM fits with mt fixed in Fig. 11(a) and 11(b) and with mt = 173.2 ± 1.34 GeV 
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(c) (d) 

Figure 11. CMSSM parameter distributions in {M^^^Myi) and (tan/3,Ao). (a) and (b) show the 
fit results for ~ (126 ± 2 ± 3) GeV and fixed mt = 173.2 GeV, compared to the LHC fit. (c) and 
(d) show the fit with the same input observable set, but with mt ~ 173.2 ± 1.34 GeV floating free 
in the fit. in comparison to the fit with fixed mt- 



floating in Fig. 11(c) and 11(d). The possible signal at rrih ~ 126 GeV shifts the allowed 
region strongly into regions of higher Mq and as compared to the LHC fit. This is 

due to the larger squark masses necessary to lift m/^ so strongly above the tree level bound 
of m/i < mz- Also, large tan/3 is clearly preferred, while again showing a flat profile in 
tan /3 outside the Icr region. As expected from the small effect which floating mt had on 
the fit in Fig. 10, there is no significant difference between the allowed parameter ranges 
for mt fixed and mt floating, albeit there is a significant jump in the best fit point due to 
the flatness of the profile. Since mt is expected to have the strongest direct effect on 
the prediction of all SM parameters, this confirms that with the current observable set the 
SM parameters can be fixed in the fit, since their uncertainties decouple completely from 
the SUSY parameter uncertainties. 

Figure 12 shows the allowed parameter space of the NUHMl model. Negative values for 
can be considered because the relevant parameter combination for EWSB is |//p + M^, 
which we checked to be positive above the electroweak scale. Since this model reaches lower 
for the same observable set as the LHC+m/j =126 fit, its area of low is deep enough 
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Figure 12. NUHMl parameter distributions in (Mp, M12) in (a), (tan/3, Aq) in (b) and (tan Mfj) 
in (c) for the fit to the full observable set including rrih — (126 ± 2 ± 3) GcV. 

with respect to the surrounding flat profile (see Section 4.1 and 4.2) to exclude the focus 
point region from the 2-dimensional 2a region. Other than that, no significant deviations 
from the CMSSM can be observed for Mq and However, the NUHMl may also 

be sensitive to exclusions from the search for — )• rr at the LHC. Unfortunately, these 
searches are not yet implemented in HiggsBounds in a model independent way, and hence 
not available for this fit. 

In Fig. 13(a), the predicted sparticle and Higgs boson mass range for the CMSSM with 
ruh = (126 lb 2 lb 3) GeV is shown. As discussed above, the high value of nih can best be 
reached for high sparticle masses, shifting up the expected squark mass scale significantly 
with respect to Fig. 9(b). The NUHMl fit is not pushed into the same area of a flat 
profile for heavy Mq and M1/2 as the CMSSM LHC+m/i =126 fit. Therefore, and because it 
can decouple the Higgs sector via the separate parameter Mfj, it prefers significantly lower 
sparticle and heavy Higgs boson mass scales than the CMSSM with the same observables. 

The observed fit results can be used to look at the range of allowed Higgs branching 
fractions in the CMSSM. These measurements are not used in the fit, and the range of 
predicted branching fraction values can be plotted within the allowed parameter range. 
Since the program HDECAY [87] is not included in the fitting process, it is used after the 
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Figure 13. Predicted distribution of sparticle and Higgs boson masses from the CMSSM fit with 
m.h = (126 ± 2 ± 3) GeV in (a) and the NUHMl fit to the same observable set in (b). 
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Figure 14. Predicted 2a range of Higgs branching fractions for the LHC fit of the CMSSM in (a) 
and the CMSSM fit with mh = (126 ± 2 ± 3) GeV in (b). The resuhs do include the theoretical 
uncertainty of the Higgs boson mass of ±3 GeV. Ratios of the potentially experimentally accessible 
branching fractions at LHC are given in in (c) and (d) for the same fits as above. 
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Figure 15. Solid lines show the results of the LHC fit (including direct detection constraints 
from XenonlOO [39]) on the {Mq, M1/2) (a) and (Ao,tan/?) plane (b). Increasing the runtime of 
XenonlOO (as aimed at with XENONlOOGoal) does not have a significant impact on these contour 
plots. Even the considerably stronger constraints that would result from increasing the target mass 
(as planned with XENONIT, shown with dashed lines) hardly affects the 95% CL regions - at the 
expense, however, of making the overall fit quality worse; the 68% CL regions, on the other hand, 
become larger in that case and the best fit point is shifted to a higher scale. 

determination of the 2a allowed parameter range to predict the ratios of lightest Higgs 
branching fractions between the CMSSM and the SM. In order to do this, the parameters 
are scanned in a 4-dimensional grid around the best fit point within their 1-dimensional 
2a uncertainties. We assume that uih can be measured with an experimental precision 
of better than 2 GeV, once the branching fractions or ratios of branching fractions are 
measured with sufficient precision. The theoretical uncertainty of the prediction for rrih 
must be taken into account for a prediction of -Scmssm/'Ssm- In Fig- 14(a) the allowed range 
of -Scmssm/'Bsm for various branching fractions is shown for the LHC fit at nih = 115 GeV. 
Fig. 14(b) shows the same for the CMSSM LHC+m/j =126 fit. In both cases one observes 
an enhancement of the bb final state. t~^t~ decreases due to different contributions to the 
tan/3 enhanced terms which need to be resummed [135, 136]. The main difference is in 
particular due to the gluino-sbottom and chargino-stop contributions in case of the bb final 
state. A significant sensitivity beyond the SM values can be observed in both cases, making 
potential measurements of the branching fractions [137] very attractive to determine the 
model parameters further and to discover a deviation from the SM even for SUSY mass 
scales beyond the LHC reach at y/s = 7 or 8 TeV. In Fig. 14(c) and 14(d) it can be seen that 
even the (at the LHC) experimentally potentially accessible ratios of branching fractions 
still exhibit significant potential for discovering differences from the expected behavior of 
a SM Higgs boson. The effect of such measurements on the fit can be expected to be even 
stronger assuming the foreseen precision at an e~^e~ linear collider (see e.g. [138]). 
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Figure 16. Predicted distribution of spartiele and Higgs boson masses from (a) the LHC and (b) 
the LHC+XENONIT fit. The fuU uncertainty band gives the 1-dimensional 2a uncertainty of each 
mass defined by the region Ax^ < 4 after profiUng over ah hidden dimensions. Note the different 
scales on the ordinate (mass) axis compared to Fig. 9. 

4.4 Implications for dark matter searches 

In the fohowing we discuss the impact of the astrophysical input parameters on our fit 
resuhs. Recah that adding the claimed signals in direct detection from the CoGeNT ex- 
periment [41, 108] to our fit leads to an unacceptably high ^X%r:cct contribution. We thus 
confirm earlier observations [139] revealing that the large scattering rates associated with 
these signals are incompatible with originating solely from neutralino DM in the CMSSM. 
The same holds for the DAMA/LIBRA [107] signal. 

Upper limits from the XenonlOO [39] experiment presently do not constrain the CMSSM 
parameters beyond the LEO, LHC (5 fb~^) data and the relic density of CDM alone. Con- 
cerning future prospects [110] for direct detection, we find that stronger constraints on 
the spin-independent scattering cross-section cjsi per nucleon that would result from an 
increased runtime (as planned with XENONlOOGoal) and a non-observation do not have 
a large impact, either. Increasing the target mass (XENONIT) and a non-observation, 
on the other hand, would enlarge the allowed 68% CL region in both the (Mo,Mi/2) and 
the (Ao,tan/3) plane (see Fig. 15) and move the best fit point to somewhat higher scales. 
At 95% CL, the contours comprise only slightly larger regions ~ with the largest change 
visible in the focus point region at Mq > 4.5 TeV. The overall fit quality, however, be- 
comes considerably worse in this case and the minimal /ndf increases from 13.1/9 to 
15.0/9, as already stated in Table 4. In fact, Fig. 16 shows that implementing the projected 
XENONIT limits has the effect of considerably raising the best-fit mass of essentially all 
particles with respect to the current baseline expectation, in particular for the heavy Higgs 
bosons and squarks. 

We also considered the combined impact of direct detection limits and a Higgs boson 
mass of rrih ~ 126 GeV. In Fig. 17(a), we only use limits from HiggsBounds to constrain 
the Higgs boson mass: in an overlay plot, we show how constraints from direct detection 
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Figure 17. Current and future limits from the Xenon experiment, when (a) using HiggsBounds 
or (b) assuming a Higgs boson mass of m^, = [126 ± 2 ± 3] GeV. In (a), the minima for the first two 
fits are identical. See the text for further comments. 



experiments will hardly improve with increased runtime (XENONlOOGoal) but would do 
so with an increased target mass (XENONIT) ~ however at the cost of an increased Xmin' 
see the comment above. Note that our XENONIT allowed region extends into the regions 
nominally excluded by the experiment. The reason for this is that we adopted the rather 
conservative choice of assigning a theoretical uncertainty of 50% to the calculation of the 
spin-independent scattering cross-section per nucleon asi, see Section 3.2.2. Restricting the 
Higgs boson mass to m/j = [126 it 2 it 3] GeV instead, we see in Fig. 17(b) that the preferred 
neutralino mass moves from 270 GeV to 497 GeV; this trend to higher masses is of course 
expected because a large Higgs boson mass generally requires rather high values for the 
SUSY breaking scale (at least in the minimal SUSY version considered here). Concerning 
future prospects for direct detection, we can see that a large Higgs boson mass worsens the 
situation as it pushes the best-fit crsi to lower values. 

We find that our indirect detection upper limits from dwarf spheroidal galaxies, us- 
ing gamma-ray observations by the Fermi satellite, are still too weak to give a noticeable 
Ax^-contribution for neutralino DM in the CMSSM. This is not a great surprise as the 
limits barely touch the annihilation cross-section of ~ 3 x 10~^^cm^/s, which is naively 
expected for thermally produced DM. Concretely, we used the photon flux upper limits 
from Ref. [140], for > 100 MeV, on neutralino pair annihilation into hh final states 
- which very often gives the dominant contribution to the total flux (mostly from pho- 
tons with E,^ <C 'm^. An improved treatment would also take into account the photons 
from other final states. However we caution that this is not straightforward to imple- 
ment in those regions of the parameter space where the photon spectra are very model- 
dependent [143]. Further improvement is possible by using updated limits from a combined 
('stacked') analysis of all dwarf spheroidal galaxy data taken by Fermi [141]. We thus ex- 
pect that a more accurate treatment of the combined gamma-ray limits in AstroFit, which 
is planned for future fits, would actually impact the CMSSM parameter space (as also found 
in e.g. Ref. [144]). This expectation is reflected in Fig. 18(a), where we show ((Jv)^^_^i^ 
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Figure 18. la and 2cr contours in the (erf vs. m-^^ plane (a), relevant for indirect dark matter 

detection, slightly change when applying current and projected limits from direct searches for dark 
matter. Also shown are the gamma-ray limits which we adopted here [140] (dashed horizontal 
line) as well as the currently most stringent limits [141] (solid horizontal line) that will be used in 
an update of this study. In (b) the annihilation cross-section vs. the spin-independent scattering 
cross-section is shown in order to demonstrate the complementarity [142] of direct and indirect dark 
matter searches. 

vs. m^j^: While the limits that we have implemented indeed do not touch the 2a regions, 
the improved limits from the joined dwarf spheroidal galaxy analysis [141] do. Those limits 
were not available in AstroFit when the scans were set up. Let us also stress that we plot 
here only the annihilation cross-section into bb final states. Future prospects for indirect 
dark matter detection are thus actually much better than what is naively inferred from 
this figure - especially when explicitly taking into account gamma-ray spectral features in 
the analysis rather than only counting the number of photons [145] . 

In Fig. 18(b), we plot the neutralino annihilation cross-section against the spin-in- 
dependent scattering cross-section, demonstrating that indirect and direct dark matter 
searches indeed probe the parameter space from an orthogonal direction [142] and are 
highly complementary even for very constrained scenarios like the CMSSM. In particular, 
improving current gamma-ray limits by about one order of magnitude (as might be rather 
straight-forward with future air Cerenkov telescopes [142]) would allow to probe models 
that are completely out of reach even for XENONIT. Models in the upper right corner of 
Fig. 18(b), on the other hand, would in principle allow for a future simultaneous detection of 
dark matter with both direct and indirect methods which evidently would make any claim 
for a corresponding signal much more convincing. We checked that adding the Higgs-mass 
constraint ruh = [126 ± 2 it 3] GeV does not have a major impact on the 2a region in this 
plane. The la region, on the other, hand blows up considerably. The best fit point moves 
to cJsi ~ 10~^^pb and {av)^,^ ~ 10~^^cm^s~^. This again just reflects the overall worse 
quality of the fit. 

The relic density of cold dark matter remains a strong constraint on the fit. Indeed, it is 
well known that only relatively small regions in the full parameter space of the CMSSM can 
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Figure 19. Solid lines show the standard LHC fit, dashed lines the impact of the relic density 
constraint on (a) the {Mq, Mi^2) plane and (b) ctsi vs. m^o. In (b) also the current (XenonlOO) as 
well as near (XENONlOOGoal) and far (XENONIT) future reaches of direct detection experiments 
arc shown (from top to bottom). 



account for thermally produced neutralino dark matter with the correct relic density: (i) 
the bulk region at low values of Mq and M1/2, which is now essentially excluded by collider 
data, (ii) the funnel region at intermediate values of Mq and M1/2 with tua — 2m^ so 
that annihilation via the pseudo-scalar Higgs boson becomes resonant, (iii) the focus point 
region at Mq ^ Mi^2 where the lightest neutralino has a sizable Higgsino component and 
(iv) the stau co-annihilation region at large but small Mq with ~ mf. Moreover, 
for Aq ^ there is the stop co-annihilation region at somewhat larger values of Mq and 
correspondingly very large values of Mi/2- The impact of the relic density constraint can 
best be demonstrated in the {Mq, M1/2) plane, as shown in Fig. 19(a). Note that for fixed 
values of Aq and tan /3 the allowed regions would only appear as thin strips in this plane. 
While we find that the rehc density calculated by DarkSUSY [86] and micrOMEGAs [82] can 
differ up to a factor two for some parameter combinations, this does not have a significant 
impact on our best fit regions and 2D contour plots. The most striking difference appears in 
Fig. 19(a) where the co-annihilation region does not extend all the way out to Mq ~ 2.2 TeV 
but only to Mq ~ 1.8 TeV when using DarkSUSY instead of micrOMEGAs. Figure 19(b) shows 
that the relic density together with the particle physics input favors neutralino masses in 
the range between 200 GeV and 500 GeV. At 95% CL an additional range between 50 and 
70 GeV is allowed. This corresponds to a resonance in the focus point region, see below. 
However, the expected spin-independent scattering cross-section is hardly affected by the 
relic density constraint. 

Note that the 2a region with m o ~ 50 — 70 GeV that appears in several of the 2D 
plots discussed above does not show up in the mass distribution plots (Fig. 16) because it 
is not contained in the 1-dimensional 2a environment. These neutralinos lie in the focus 
point region appearing in Fig. 19(a) at very low values of M1/2 and Mq ~ 4 TeV. They have 
a Higgsino fraction of roughly 1%. Their relic density is almost exclusively set through 
s-channel annihilation via the light scalar Higgs h. 
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Moving along the focus point 'strip' to larger masses, the Higgsino fraction increases 
and at M;^/2 ^ 2.5 TeV the neutrahno becomes an essentially pure Higgsino that acquires 
the correct relic density for m^o ~ 1.1 TeV. In our fit, however, models with such large 
neutralino masses do not appear because the sampling algorithm runs out of statistics 
at these very large values of Mq > 4.5 TeV due to the excessive inter-parameter fine- 
tuning (discussed in Section 4.6) for such high Mq^. For the same reason, we do not 
sample the tail of the co-annihilation strip which in principle admits TeV-scale neutralinos 
co-annihilating with almost degenerate staus. We mention that the very small region at 
Mq ~ 4 TeV and slightly higher values of M1/2 also belongs to the focus point region. Albeit 
in this case it contains models where the lightest neutralino has a sizable Higgsino fraction 
of 0(10%) and a mass around ~200GeV. The relic density is thus mostly set through 
(non-resonant) annihilation into SU{2) gauge bosons, which is helped by the presence of 
relatively light charginos and neutralinos, though not light enough that co-annihilations 
would be important. 

4.5 Fine-tuning 

A major motivation for supersymmetry at or near the electroweak scale is to solve the 
hierarchy problem [147, 148]. This is a fine-tuning problem in the Higgs sector of the 
theory, where the low-energy observable Higgs boson mass depends sensitively on the input 
parameters at a postulated new scale of physics, such as the unification scale Mqut? or the 
Planck scale. This fine-tuning is removed and the hierarchy problem solved for unbroken 
supersymmetry. A solution is retained after supersymmetry breaking, if the supersymmetry 
masses are not too heavy. How heavy depends on how much fine-tuning is perceived to be 
acceptable. This is the origin of the expectation, that the super symmetric masses should 
be MsusY < 0(10 TeV), or even less. 

Beyond the Higgs boson mass, also other low-energy parameters can be fine-tuned, 
such as the electroweak breaking scale, represented by the boson mass, mz- The fine- 
tuning can be quantified by various measures. A popular one was introduced in [79]. It is 
based on the logarithmic derivatives of any given observables with respect to the parameters 
considered. For mz considered as a function of variables Oj the relevant numbers are 

Oj dm\{ai) 



ciai = 



m 



(4.1) 



The overall fine-tuning at a given parameter point is then given by 

A = maxj [c(aj)] . (4.2) 

It has since been widely used to quantify the fine-tuning in the (C)MSSM both before 
[149-154] and more recently also including LHC data [155-157]. 

For the CMSSM, SOFTSUSY [81] can be configured to calculate ah of the derivatives in 

1/2 



Eq. (4.1). The set of parameters {ui} consists of Mq, M1/2, Aq, ms, /x, and ht-^ Points 



* See, however, Ref. [146] for a detailed discussion of the focus point region arguing that, even for very 
large values of Mq, the actual fine-tuning is small if expressed in physical parameters. 

^SOFTSUSY also allows the computation of the derivatives with respect to the non-CMSSM parameters 
ma , the soft breaking gluino mass, and ht , the top Yukawa coupling. 
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Figure 20. The minimal amount of fine-tuning as a function of Mq and M1/2, where Aq and tan/3 
are profiled, for (a) the pre-LHC fit and (b) the LHC fit. (c) shows the difference, (b)— (a), between 
both fits, which is mostly close to vanishing in the overlap of the 95%-CL regions. The last bin 
on the z-axis in (a), (b), and (c) contains all points with A > 500. The smallest values for the A 
parameter are 39.7 (193.5) for the pre-LHC (LHC) fit. 

with larger values for A are more fine-tuned and are in general regarded as less natural. A 
fit has been performed to assess the fine-tuning, using the same observables and LHC limit 
as for the LHC fit before, but with SoftSUSY as spectrum calculator instead of SPheno. 

As we have seen, the LHC exclusion leads to an overall decrease of the fit quality. 
Despite having a more constrained system, this worse fit can lead to more points in the 
2a allowed range around the best-fit point. In particular, for a fixed value in the Mq-Mi/2 
plane a larger range of Aq and tan /? values can be allowed. In Fig. 20 we plot the lowest 
values of A in the Mq-AIi/2 plane, where the parameters Aq and tan/3 have been profiled. 
In Fig. 20(a), A is shown for the pre-LHC fit. The lowest amounts of fine-tuning (> 39.7) 
are observed for low values of Mq and M1/2, while in particular for high values of M1/2 the 
A-parameter reaches larger values, over 300.^ 

Fig. 20(b) shows the fine-tuning taking into account the LHC exclusion with a dataset 

®We also considered the fine-tuning as a function of Mq and M1/2 for the best-fitting values of Ao and 
tan /3 in each bin. The results are indistinguishable from Fig. 20(a) . 
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of 5fb~ . This fit introduces points witfi a significantly larger fine-tuning (A > 500) 
compared to the pre-LHC fit and the point of lowest fine-tuning (A = 193.5), is significantly 
more fine-tuned than in the pre-LHC fit. This is because the fit pushes one to higher values 
of Mq, Ml 1 2- Again overall in the LHC fit, low values of Mq and allow for a lower 

fine-tuning, while for larger values for both parameters, in particular in the focus point 
region, a higher fine-tuning is observed (see also discussion in Section 4.4). 

Fig. 20(c) shows the difference, (b) — (a), between the two fits in the region that is 
included in the 2o"-contours of both the pre-LHC as well as the LHC fit. The solid/dashed 
lines indicate the part of the parameter space which is not excluded at the 95%-CL in 
the pre-LHC/LHC fit, but which is excluded at the 95%-CL in the LHC/pre-LHC fit. In 
the (Mo,Mi/2) plane, only a few points with a significant change in the fine-tuning show 
up, while for the majority of points which are not excluded in both fits, the fine-tuning is 
comparable. 

We find slightly larger minimal values of A compared to other studies, for instance 
[155]. However, the values we find are not dramatically larger and our fits are constrained 
by a larger set of observables, such that deviations are expected. While we find that with 
earlier LHC limits (£int ~ lfb~^), using the criteria of [155], points with an acceptable 
amount of fine-tuning (A < 100) remain in reasonable agreement with all precision mea- 
surements, the minimal value of A is pushed above this threshold by the up to date LHC 
results. 

4.6 Correlation of fundamental parameters 

The fine tuning measure A described in the previous section has many merits from a 
theoretical point of view. However, from an experimental point of view, there is one major 
shortcoming: The absolute change of an arbitrarily selected observable such as mz is 
experimentally of limited interest. The experimentally more meaningful quantity is the 
change of the observable relative to its uncertainty. Furthermore, experimentally we are 
interested in all observables chosen to calculate the x^, not only in one particular quantity 
of special theoretical interest. 

Therefore, we propose a different view on looking at the naturalness of a model in 
the light of all observables and all uncertainties, inspired by the fine tuning. There, if one 
prediction for an observable depends strongly on varying one parameter, the agreement of 
all measurements with the predictions must be restored by varying other parameters in a 
correlated way. Therefore, we examine the correlation of the fundamental parameters in 
the fit as an additional information about how highly the parameters have to be tuned 
with respect to each other, measured in the form of the quantity ^max defined below, in 
order to fulfill the same level of agreement with all data in the fit. Thus, all observables in 
the fit contribute with their current experimental precision to fPmax- 

A priori, the fundamental parameters of the tested models are independent. If however, 
starting from a point in parameter space with a quality of fit x^, if one parameter is varied, 
the experimental bounds may force a correlated change in one of the other parameters, 
in order to maintain an equally good fit. A first order measure for the dependence is 
Pearson's product-moment correlation coefficient. For each pair of parameters, Pj, Pj € 



- 35 - 



{Mo, Ml/2, ^Oi tan/3}, it is defined by 

\ (^P.<^P, I 

{Pi) is the average value for the parameter Pi and a p. the standard deviation of the Pa- 
rameter Pi. Therefore \Qij\ < 1. 

The set of values of the parameters {Pi} used for the computation of the average and 
standard deviation at a fixed point P in the parameter space is defined by the value of 
at this point, Xp- We consider all points in the MCMC scan with a quality of fit in the 
range [xpjXp + where Ax^ = 0.001. These are very thin slices. All points included 

twice or more often in the MCMC scan are taken into account only once. If less than 6 
points exist within the Ax^ < 0.001 boundary the tested point is assigned no value and 
excluded from the further analysis. 

For each point in the 95%-CL region, we now compute the quantity 

^•max = maXjj {\Qij\) , (4.4) 

considering all 6 combinations of {Pi, Pj). For a fixed value of Mq and M1/2, and with Aq, 
tan/3 profiled, we then compute the minimum of ^max, [f?max]min, over the 95%-CL region. 
Large values of [^max]min mean the fit is highly correlated. 

A model highly constrained by the observables with their given precision will nat- 
urally yield ^max ^ 1 near the point with the exact minimum. Depending on the 
model, observables and precision, it may however show high {e.g. near 1) or low values 
of [£'max]min over the full allowed range, depending on how much the parameters need 
to be fine-tuned with respect to each other in order to achieve a given agreement between 
the model predictions and all data, taking the experimental and theoretical precision fully 
into account. 

[^max]min Can describe features of the analyzed parameter space at very different scales. 
For low values of its value is dominantly built from the small scale structure of the 
parameter space. For larger values of x^, all parameters vary over a larger range and thus 
the included parameter space grows significantly and correlations on a more global scale 
are reflected. 

In order to have a baseline comparison, we have first implemented this procedure 
without the experimental constraints considered throughout this paper. Instead we have 
explored the correlation where only a minimum set of requirements has been used to 
constrain the parameter space: we require (a) the lightest neutralino to be the LSP, (b) 
consistent radiative electroweak symmetry breaking and (c) the absence of tachyons. In 
order to have a meaningful comparison, the correlation between two parameters is not 
calculated over an infinite parameter range, but instead in a rectangular subset of the 
parameter space, the boundaries of which are deflned by the maximum and minimum 
values of each of the four parameters as found in the 95%-CL region of the LHC fit of the 
same model. After performing this procedure we find the maximum correlation between 
parameters with these minimal constraints is found to be smaller than 0.2. This is what 
the following results should be compared against. 
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Figure 21. Minimum value of gmax as a function of Mg and A/i/2> where and tan/3 are profiled, 
(a) shows the minimum correlation in the 95%-CL region of the pre-LHC fit, (b) shows the minimum 
correlation in the 95%-CL region of the LHC fit, (c) shows the difference of both distributions. For 
nearly all values of Mq and M1/2, the pre-LHC fit shows a larger minimum correlation. The smallest 
values of g^ax are 0.62 (0.29) for the pre-LHC (LHC) fit. 



The results of this study including the pre-LHC and LHC constraints are summarized 
in Fig. 21. In Fig. 21(a) we show [iPmaxlmin in the M0-M1/2 plane, when including only the 
pre-LHC constraints. ^0 and tan/3 are profiled. We see that due to the now restricted 
parameter space the correlations are all above 0.75, and well above the 0.2 considered 
above. Moreover the correlation increases when increasing Mi 12 and in particular when 
increasing Mq 

In Fig. 21(b) the LHC limits are now also taken into account. The colored region 
is significantly extended because the quality of the fit has gone down and the 2cj region 
is larger. At the same time, because ^0 and tan /3 can extend over wider ranges, the 
correlation is markedly decreased compared to Fig. 21(a). It is however still everywhere 
above the baseline value of 0.2. 

In Fig. 22, we consider slightly modified scenarios, with one additional fit parameter. 
In all cases we have fixed the lightest Higgs mass to rrih = 126 GeV. In Fig. 22(a) [Fig. 22(b)] 
we show the minimal correlation for a fixed (free) top mass, m^. Due to the fixed Higgs 
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Figure 22. Minimum value of gmax as a function of Mq and where and tan /? are profiled, 

(a) and (b) show the minimum correlation for a fit of the CMSSM model with a fixed (free) mj. 
For all points, the minimum correlation increases for a free nit. The minimum value for [fJmaxlmin 
is 0.59 (0.70) for the fit with fixed (free) m*. In addition, (c) shows the ft„ax as a function of Mq 
and My 11 for a fit of the NUHMl model, where Aq, tan^ and TO//o are profiled. The minimum 
value of gimax for this fit is 0.45. All fits use mh ~ 126 GeV. 

mass the quality of the fit degrades significantly. The 2a region is thus now extends to very 
high values of Mq. For the free top mass the correlations are slightly higher than for the 
fixed case. In Fig. 22(c) we show the correlations in the case of the previously described 
model NUHMl. The minimum value of the correlation for this fit is 0.45. In all cases 
shown in Fig. 22, Aq and tan /3, as well as nit and M/, where applicable, are profiled. 

As discussed above, fPmax may vary significantly for similar values of Mq and if 
additional points with worse agreement between model predictions and measurements are 
reintroduced by adding another observable. This was so far taken into account by letting 
vary over the full accessible range, and then computing [f?max]min- It is also interesting 
to simply consider the best-fit value of for fixed values of Mq, Mi For the fits of 
the CMSSM model with rrih = 126 GeV, this is shown in Fig. 23. Figs. 23(a)/23(b) show 
Qma.x in the M0-M1/2 plane at the point with the lowest x^ P^r bin for the mf=fixed/mt 
= free fit. Figure 23(c) shows the difference between the two fits. A negative (positive) 
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Figure 23. The value of Qmax as a function of Mq and M1/2 at the point with the lowest x^, where 
^0 and tan/3 arc profiled, (a) shows the minimum value for a fit of the CMSSM model where rrit 
is fixed, (b) shows the minimum value for a fit of the CMSSM model where mt is free, (c) shows 
the difference of the two distributions. Near the best-fit point, the minimum of ^Pmax is increased if 
rrit is fixed, while at the boundaries of the 95%-CL region the minimum value is decreased. Both 
fits are constrained by mh = 126 GeV. 

entry indicates that with mt=free(mt=fixed) a lower correlation between the parameters 
is found. In particular close to the minimum, nit = free allows for a lower correlation, 
while in other regions the maximum correlation between the varied parameters increases if 
the top mass is free. White areas indicate that at the corresponding point with the lowest 
there were less than 6 points within Ax^ < 0.001 such that no meaningful correlation 
could be calculated. 

4.7 Including Q as a nuisance parameter 

When fitting constrained models such as the CMSSM, the SUSY parameters are fitted 
at the GUT scale. The predictions for the observables are calculated from the MSSM 
Lagrangian calculated at a scale of Q = 1 TeV where we fixed the scale according to 
the SPA conventions [158]. Theoretical uncertainties on the predictions are included in the 
calculation of the &s noted in Section 3. However, this does not included the parametric 
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Figure 24. Fits including LHC exclusions, with and without Q as a nuisance parameter. Details 
about the best fit point are given in Tab. 4. Results are shown in the (A/q, M1/2) plain. The upper 
limit on the SUSY parameter scale is strongly lifted upwards once the RGE uncertainty is explicitly 
included in the fit by floating the scale Q. 

uncertainties arising from the fact that the masses are calculated using two-loop RGEs and 
one-loop (two-loop) threshold corrections for the sparticle (Higgs) masses. It is generally 
assumed in previous work in the literature (see e.g. [64]), that these uncertainties are 
negligible or at least sub-dominant for the observables described in Section 3. 

However, there is a straight-forward possibility to include these previously neglected 
uncertainties: The scale Q, where the spectrum is calculated, can be included as a nuisance 
parameter into the fit. It is allowed to float within ^ yj'^^li ''^h < Q < 
the parameter uncertainties on the physical parameters of the model can be obtained by 
profiling over the hidden dimensions including Q. As shown in Fig. 24 for a fit including 
LHC exclusion, HiggsBounds limits and the full AstroFit limits, this has a sizable effect 
on the upper limit of the available parameter space, which is dominantly determined by 
{g — 2)^ and ^dm- The reason for this large increase is the flatness of the distribution 
which implies that already a small change Ax^ can lead to the observed changes for the 
allowed regions. However, the best regions hardly changes. This result suggests that in 
future fits, the RGE uncertainty should be included as well. 

We note for completeness that for very large values of Mq and/or M1/2 one obtains 
1 TeV < ^ -^"^ti "^42 ■ However, in these regions the theoretical uncertainties for the mass 
calculations in codes likes SPheno or Sof tSUSY is most likely larger than the scale variations 
suggests because here one should use a multi-scale approach as discussed in [159-161] in 
contrast to the single scale approach of the current implementations [80, 81]. 

5 Conclusions 

We have presented a global fit of the constrained MSSM and a non-universal Higgs mass 
model (NUHMl) to low-energy precision measurements, dark matter observables and, in 
particular, current LHC exclusion limits from direct SUSY searches in the zero-lepton plus 
jets and missing transverse energy channel. A restricted model with universal soft SUSY 
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breaking parameters at the unification scale like the CMSSM does not provide a good 
description of all observables. The best fit predicts typical sparticle masses in the TeV- 
range, albeit with large errors. We find that the description of the low-energy observables, 
{g — 2)^ in particular, and the non-observation of SUSY at the LHC become increasingly 
incompatible within the CMSSM. Note that the LHC exclusion in the zero-lepton, jets 
plus i?™'^^ channel mainly constrain the squark and gluino masses, while {g — 2)^ and 
other low-energy observables are sensitive to the color-neutral sparticles. Super symmetric 
models with common scalar and gaugino masses like the CMSSM connect these two sectors 
and thus lead to an increasing tension between {g — 2)^ and the LHC search results, which 
would not necessarily be present in more general SUSY scenarios. 

In general, there is an intricate interplay between the different contributions to the 
global fit, and an accurate estimate of the LHC exclusions also in parameter regions away 
from the published 95% CL contours is mandatory. We have thus employed an elaborate 
simulation of recent LHC search analyses based on state-of-the-art Monte-Carlo tools and 
a public detector simulation. We note that it appears challenging to provide a fast and 
accurate implementation of LHC results in a global fit for more general SUSY models with 
a larger number of parameters. 

A very strong constraint on SUSY can be expected from a measurement of the light 
Higgs boson mass. Within the global CMSSM fit the light Higgs scalar is Standard- 
Model-like, with a preferred mass around nih ~ 117 GeV. We have studied the impact of a 
potential Higgs observation on the global fit and find that a light Higgs boson with mass 
mil ^ 125 GeV is hardly compatible with the other data within the CMSSM. The large 
sparticle mass scales needed to accommodate Higgs bosons with m/j > 125 GeV drive the fit 
into a region which is incompatible with {g — 2)^ and i?-physics observables. This tension is 
reduced in models like the NUHMl where the Higgs sector is decoupled from the squark and 
slepton sector. Still, even in the NUHMl, the interplay of B{Bs — )• /u/i) and the lightest 
Higgs boson mass provides strong constraints. We have argued that important further 
constraints on SUSY models can be expected from improved limits on B{Bs — )• ^^u) and, in 
particular, a measurement of the Higgs branching fractions or ratios of branching fractions. 
For the NUHMl, the search for pseudoscalar Higgs bosons in the channel ^4*^ — )• rr will be 
another interesting input. 

We have discussed in detail the implications of the CMSSM fit for direct and indi- 
rect dark matter searches. While direct dark matter detection claims from CoGeNT or 
DAMA/LIBRA cannot be accommodated within the CMSSM, the current upper direct 
detection limits provided by XENON do not restrict the CMSSM parameter space beyond 
the constraints from low-energy observables, LHC limits and the dark matter relic den- 
sity. Projections of the direct dark matter detection limits based on a non-observation by 
XENONIT, on the other hand, would provide a further important probe of the available 
SUSY parameter space. The current indirect detection upper limits as implemented in 
AstroFit are still too weak to have an impact on the CMSSM fit. However, a more com- 
prehensive future treatment of indirect detection limits may have the potential to probe the 
SUSY parameter space in regions which are complementary to direct detection constraints. 

In addition, we have studied the fine-tuning within the SUSY models in two different 
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ways. For the fine-tuning parameter A, measuring the relative dependence of mz on all 
parameters, we find a strong increase between the fit without and with LHC direct search 
limits. For the CMSSM, a minimal value of A = 193.5 is found for the LHC fit, which is 
generally regarded as unattractive. In addition, we introduce a new concept in the form 
of the relative precision at which the parameters need to be fine-tuned with respect to 
each other, in order to achieve the same overall agreement with all relevant observables 
relative to their uncertainties. Here, we find that within the 2(t uncertainty range the 
maximal correlations amongst parameters are still in the range of Pmax ^ 0.6, hinting at 
an acceptable level of correlations. However, the correlation naturally increases strongly 
at the very minimum of the fit. 

To further explore the implications of the LHC searches for supersymmetry it is manda- 
tory to include a larger set of LHC signatures, including those which directly probe the 
non-colored sparticles, and to consider global fits of less restricted SUSY models beyond 
the CMSSM. The challenge for such future global analyses will be to provide a fast and 
accurate implementation of the LHC results, if possible for an accurate combination of 
diff'erent search channels, for more general SUSY models. 
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